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Einstein's theory of general relativity is written in terms of the variables obtained from a 
conformal-traceless decomposition of the spatial metric and extrinsic curvature. The determinant of 
the conformal metric is not restricted, so the action functional and equations of motion are invariant 
under conformal transformations. With this approach the conformal-traceless variables remain free 
of density weights. The conformal invariance of the equations of motion can be broken by imposing 
an evolution equation for the determinant of the conformal metric g. Two conditions are consid- 
ered, one in which g is constant in time and one in which g is constant along the unit normal to the 
spacelike hyp ersurf aces. This approach is used to write the Baumgarte-Shapiro-Shibata-Nakamura 
system of evolution equations in conformally invariant form. The presentation includes a discussion 
of the conformal thin sandwich construction of gravitational initial data, and the conformal flatness 
condition as an approximation to the evolution equations. 



I. INTRODUCTION 

The conformal-traceless decomposition of the gravitational field was originally introduced by Lichnerowicz and 
York 0, 13 in their work on the initial value problem. Since that time the same decomposition has appeared in various 
guises in mathematical and numerical relativity. The BSSN 0, (Baumgarte-Shapiro-Shibata-Nakamura) system 
of evolution equations for general relativity is currently in widespread use in the numerical relativity community. 
The BSSN equations are based on a conformal-traceless splitting. The recently discovered conformal thin sandwich 
construction 1^ Q , like the earlier techniques for solving the initial data problem, uses the conformal-traceless de- 
composition. Unlike its predecessors, the conformal thin sandwich scheme mixes a subset of the evolution equations 
with the constraints to facilitate the analysis. The conformal flatness condition (CFC) is an approximation to the 
Einstein equations that uses the conformal-traceless decomposition. It has been applied to numerical simulations of 
binary neutron star systems and supernovae j^. 

The conformal-traceless (CT) decomposition is 

hab = f'^gab , Kab = <^"^ + ^'P'^gabT , (1) 

where hab is the physical spatial metric and Kab is the physical extrinsic curvature. The definitions are redundant, 
in the sense that multiple sets of fields gab, -A-ab, t yield the same physical metric and extrinsic curvature. To be 
precise, these definitions are invariant under the conformal transformation^ 

gab — > gab = S.'^gab , (2a) 
<P — ^ <^ = r V , (2b) 

Aab > Aab = C'^Aab , (2c) 

r — > f = r , (2d) 

for any field ^. In many recent applications involving the conformal-traceless splitting the conformal invariance ^ 
is broken, and the redundancy in the CT variables removed, by the condition that gab should have unit determinant. 
There are various insights to be gained by leaving the determinant of gab unspecified at this point in the analysis. 
Thus, / will not impose g = det{gab) = 1. 

In most of the original works on the initial value problem and on gravitational degrees of freedom 0, ll(Tj | the 
determinant g is not restricted. Rather, the field gab is described as representing a conformal equivalence class of 
metrics that defines a conformal geometry. The approach taken here is equivalent. However, instead of describing 
the results in terms of conformal geometry, I treat gab as an ordinary spatial metric (symmetric type (2) tensor) and 
emphasize the conformal invariance Q of the action and equations of motion. 



^ The arrow notation means the following: Replace the fields g^b, V'l ^uti and r with their "barred" counterparts, which differ from the 
unbarred fields by certain factors of ^. In Eq. Q the factors of ^ completely cancel in each term so that hab and Kab are independent 
of ^. Thus, hab and K^b are conformally invariant. 
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The physical metric and extrinsic curvature are also invariant under the "trace transformation" defined by r — > 
f = T + C and Aab — <■ Aab = Aab — Cv'^ffah/S, with gab and tp unchanged, for any field C- This redundancy in the CT 
variables is usually removed by the condition that Aab should have zero trace. It turns out that there is no particular 
advantage to be gained by leaving the trace of Aab unspecified. Thus, / will impose A = Aabfj""^ = 0. Then the CT 
variable t is the trace of the extrinsic curvature, r = Kabh"'^ ■ 

At first sight it might appear surprising that I have chosen to keep the conformal invariance but eliminate the 
trace invariance. The reason is that the trace transformation is essentially trivial. It just corresponds to a change in 
the splitting of the trace of Kab into the two terms r and Alp~^ . With the choice A = 0, the trace of the extrinsic 
curvature is placed entirely in r. In a similar fashion, the conformal transformation can be viewed as a change in 
the splitting of the determinant of hab into the two factors ip^"^ and g. However, this is a nontrivial change precisely 
because the field ^ can carry a nonzero density weight. Throughout this paper I define the "unbarred" CT variables 
that appear in the definitions as having no density weight. Thus, gab and Aab are spatial tensors of type (2), and 
Lp and r are spatial scalars. If ^ carries a density weight, then the conformally transformed variables, the "barred" 
variables in Eqs. (j^J, acquire density weights. 

By not specifying <? = 1 from the outset, the equations of motion we obtain are conformally invariant. We can 
then consider breaking the conformal invariance by fixing g. In principle g can be chosen as any t-dcpendent scalar 
density of weight 2. In practice I expect that the most natural way to break conformal invariance is to choose g 
arbitrarily on the initial time slice and then evolve g according to some simple prescription. Two natural prescriptions 
are considered, one in which g remains fixed along the time flow vector field (the "Lagrangian condition") and one in 
which g remains fixed along the normal to the spacelike hypersurfaces (the "Eulerian condition"). The equations of 
motion differ between these two cases in the way that the shift terms appear. These differences can be interpreted as 
changes in the density weights of the CT variables. 

In this paper I carefully examine the full Einstein theory, constraints and evolution equations, in terms of the 
conformal-traceless variables. One of my motivations for this analysis is to help clarify various subtle issues that 
arise along the way. Most of the subtleties concern density weights. Consider what happens when 17 = 1 is fixed from 
the outset. Then the conformal metric and conformal factor must carry nonzero density weights, because otherwise 
g would not remain equal to 1 under a change of spatial coordinates. This leads to some unusual and potentially 
confusing results. For example, in the initial value problem the Hamiltonian constraint is written in terms of the 
conformal Laplacian operator acting on the conformal factor. Yet, the conformal Laplacian acts as if neither the 
conformal metric nor the conformal factor carries any density weight. It is clear that the experts on the initial value 
problem already recognize and understand this subtle point. One of my goals for this paper is to present a framework 
for the conformal-traceless decomposition in which these subtle issues are largely avoided. 

A second motivation for this analysis is to help pave the way for the development of numerical techniques in general 
relativity based on the "variational integrator" approach to numerical modeling. (For an overview, see Ref. |Tl|.) In 
this approach, finite difference equations are derived by extremizing the action functional in discrete form. It turns 
out that for a constrained system like general relativity, the equations obtained from a discrete variational principle 
will mix the constraints and the evolution equations at each timestep. If the fundamental variables are the physical 
metric and extrinsic curvature, then one is faced with the problem of solving a discrete version of the thin sandwich 
equations jl^ . By now it is well known that the thin sandwich equations are ill-posed, and do not admit a solution 
in the generic case. On the other hand, the conformal thin sandwich equations are expressed in terms of "nice" 
elliptic operators and one expects the system to be generically well-posed. In order to apply the variational integrator 
construction to general relativity, we must first express the action in terms of the conformal-traceless variables. 

Throughout this paper I use the lapse anti-density a as the undetermined multiplier for the Hamiltonian constraint. 
The lapse anti-density carries density weight —1 and the Hamiltonian constraint carries density weight 2. York and 
collaborators |a, 0, Il3j have pointed out a number of reasons why this choice is preferred over the traditional scalar 
lapse function. It has also been shown that the BSSN equations are equivalent to a strongly hyperbolic system 
with physical characteristic speeds when the lapse has density weight —1. 

In Sec. II I begin by writing down the action and equations of motion in Hamiltonian form. In Sec. HI the action 
is expressed in terms of the CT variables Because the action is conformally invariant, it can also be expressed in 
terms of a set of conformally invariant variables that includes the density weight —2/3 metric gab ~ g^^^^9ab- The 
equations of motion for the CT variables are derived in Sec. IV. There, it is pointed out that the conformal invariance 
can be broken with either the Lagrangian or Eulerian condition on the determinant of gab- This leads to two sets of 
evolution equations that differ in the way that the shift terms enter. In Sec. V I examine the BSSN system. It is 
written first in conformally invariant form, then with the Eulerian and Lagrangian conditions applied. The conformal 
thin sandwich equations are derived in Sec. VI. They follow directly from the Hamiltonian and momentum constraints 
and the evolution equations of Sec. IV. The CFC approximation is discussed in Sec. VII. There it is pointed out that 
the CFC equations are identical to the conformal thin sandwich equations with the restrictions that the conformal 
metric gab is flat and the time slicing is maximal, r = 0. The main results are summarized briefly in Sec. VIII. 
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II. ADM ACTION AND EQUATIONS OF MOTION 

The Arnowitt-Deser-Misner (ADM) gravitational action is 

S^'^^ [hab, a, = [ dt d^x {P^^hab - aH^^^ - P^Mi^A , (3) 



which is a functional of the physical 3-metric hab, its conjugate P"'', the lapse anti-density a, and the shift vector 
The dot denotes a time derivative, hab = dhab/dt. The gravitational momentum is related to the extrinsic curvature 

by 

P^b = lVh{h''^h'"^ - h'"'h'"^)Kcd ■ (4) 

The gravitational Hamiltonian and momentum densities are defined by 

^(g) ^ 2P"''F„f, - p2 _ , (5a) 

Mi^^ EE -2VbP'a , (5b) 

where P = P'^^hab is the trace of the gravitational momentum and h is the determinant of hab- Also, Va and 
TZ = h°'^TZab arc the covariant derivative and scalar curvature for the spatial metric hab- If matter fields are present 
the complete action S = 5^^^ + 5'^'"^ includes a functional S*'™^ of the matter fields in addition to the gravitational 
action. I will assume that the matter is "minimally coupled" to gravity, so the matter action does not depend on 
derivatives of the spatial metric or on P'^''. Throughout this paper I use units in which SttG = 1, where G is Newton's 
constant. 

Variations of the ADM action 10) plus matter action i'f™) with respect to the lapse anti-density and shift vector 
yield the constraints 

n^3'> = -hp , Mi^^ = Vhja , (6) 

where p and ja are the energy and momentum densities for the matter fields. They are defined, along with the spatial 
stress tensor s°^, in terms of the functional derivatives of the matter action by 

h Sa ^ ° ^//i Sf3°- ' ah Shab h 5a 

Variations of the action with respect to hab and P"'' yield the well-known ADM evolution equations 

d^hab = a(^4Pab - 2Phab) , (8a) 

a^LP"'' = ai-AP'^'Pl: + 2PP°'' + -hUh'"'' - -hn''^\ + -hX/^XI^a - -hh'^^V'^a + -ahis"^ - ph'''') . (8b) 
V 2 2/2 2 2 

The derivative operator on the left hand sides of these equations is defined by 

di_ = d/dt - Cp (9) 

where £^ is a Lie derivative along the shift vector field /3°. Thus, {a\/h)~^dj_ is the derivative with respect to proper 
time along the unit normal to the spacclike hypersurfaccs. Specifically, the Lie derivatives are defined by 

Cfshab EE P^d.hab + hacObP' + h.bdaP' = 2V(„/3fc) , (10a) 

CpP"^ EE dciP^P"'') - P^^dcP^ - P^^dcP" = Vc(/3"P°^) - 2P"(°Vc/3''^ , (10b) 

where da = d jdx^'' and (a • • • 6) denotes symmetrization on the indices a and h. Note that hab is a symmetric type (2) 
tensor on the spatial manifold, and P°'^ is a symmetric type (g) tensor density of weight 1 on the spatial manifold. 
(The terminology of density weights is fixed by noting that \fh is a scalar density of weight 1.) 

The time derivatives in the action ||2J) can be written in terms of the normal derivative 9^ ■, with the result 

S'(s)[/i„fc,P'''',a,/3"] = [ dtd^x\p''''d±hab-aH^s'>] , (11) 



The difference between the actions (O and (|11() is a total derivative, 2 J dt d^ xV aiP""^ Pb) ■ This term vanishes for 
compact spatial manifolds, and integrates to a boundary term for non-compact spatial manifolds. The presence of 
boundary terms in the action affect the admissible boundary conditions to be imposed in the variational principle, 
but they do not affect the resulting equations of motion. I will ignore such boundary terms in the present work. 
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III. ACTION IN CONFORMAL-TRACELESS VARIABLES 



With the physical spatial metric and extrinsic curvature written as in Eq. the gravitational momentum Q 
becomes 

P"' = -\^gv~'A^' + ivs^Vs"" . (12) 

Recall that Aab is trace free, ^ = 0. Also note that indices on Aab are raised with the inverse conformal metric g"'^ . 
We can express the ADM action in terms of the conformal-traceless (CT) variables by substituting the expressions 
for hah and P"'' into Eq. ((TTl)- The result is 



(13) 



where the gravitational contribution to the Hamiltonian constraint is 



^ ^gA-'Aa, ~ ig^iV^ - ig^^i? + ig^^D^^ . (14) 

Here, R is the scalar curvature of gat-, and Da is the covariant derivative compatible with gab- Note that the "velocity" 
term d± {g~^^^gab) in the action is trace free. It can be written as d±_ {g~^^^gab) = g~^^^ {d±gab)'^^ , where TF stands 
for the trace-free part of the expression enclosed in parentheses. 

Observe that the traceless property of Aab must be preserved when the action is varied. This condition can be 
enforced with a Lagrange multiplier. The result is equivalent to simply demanding that the functional derivative 
with respect to Aab must be traceless. In general, if we have a functional S'[Tat] of a traceless tensor Tab said 
SS = J d'^x E'^''STab, then extremization of «S'[Tafc] yields (£'"'') = 0. This result, that the functional derivative 
with respect to a traceless tensor is trace-free, should not be too surprising. For example, we routinely vary action 
functionals with respect to symmetric tensors. The symmetry property can be enforced with a Lagrange multiplier, 
but the result is equivalent to simply demanding that the functional derivative should be symmetric. 

A key piece of the calculation for the action (|13|l is the expression for the scalar curvature of the physical metric 
hab in terms of the conformal metric gab and conformal factor (p: 

n = ip-^R-8ip-^D'^ip . (15) 

This result can be derived in a straightforward way by inserting the decomposition hab = 'fi'^gab into the definition 
for the Ricci scalar, 7^ = 2g°'^ {d[cXl}a + ^d[c^b\^ ' ^^ere FJJ^ = h'"^ {dbhdc + dchd - ddhc) /2 are the Christoffel 
symbols. 

The gravitational field contributions to the momentum constraint are "hidden" in the d± terms in the action, since 
these terms include Lie derivatives with respect to the shift vector The Lie derivatives are defined by the tensor 
character of the variables on which they act. For example, Cpgab and CpA""^ are given by expressions H1U|I with hab 
replaced by gab, P""^ replaced by A"''^ , and Vq replaced by Da- (Indices on /?" and Dq, like Aab, are raised and lowered 
with the conformal metric gab and its inverse g'^^ ■) In the calculations here and below, it is convenient to use the 
formula 

CfiT = g^'^Cp (rg-"'/^) + wTD,f3' (16) 

for the Lie derivative, where T is a tensor density of weight w. (The indices have been suppressed on T.) Then the 
shift terms in the action H13|l are 



. dtd^x^\A''''{Daf3bf^ - IrDa (/3V')1 , (17) 

shift terms / L o J 



and the gravitational contribution to the momentum constraint is 

Mi^^^-^ = ^gDbAl^-^g^'DaT . (18) 

The full Hamiltonian and momentum constraints are = 7i'^' + gtp^'^p and = A^i^'' — y^^p^ja, respectively. They 
are invariant under the conformal transformation jSJl. 
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The action as written in Eq. H13|l is a functional of the symmetric type (2) tensor gab^ the symmetric type (2) 
tensor Aab-, the scalars if and r, and the lapse anti-density and shift vector. Because the CT variables are redundant 
the equations of motion obtained by varying the action are not independent. To be precise, the action is conformally 
invariant so its variation vanishes when the CT fields are varied by an infinitesimal conformal transformation 
This leads to the relation 



= 4gafc 



Sgab 



Sip 



' SAnb 



(19) 



for all field configurations, not just those that satisfy the classical equations of motion. 

The conformal invariance of the theory defined by the action H13|) can be displayed in an elegant form by treating 
it as a constrained Hamiltonian system [l^. In that case the action becomes 



Sis) 



dt d-^x 



' metric ^- 



Q metric 
ab 



Pphi 9± QP'" + T'rootg 3^2™°*^ - a7^(^) - eC 



(20) 



where 'P^otric ^phh and 'Prootg arc the momenta conjugate to the canonical coordinates QJJ^°*''"^ = g^^^^gab, Q^^' — ^ 
and = respectively. These variables are restricted by the first class constraint C = QP^'Pphi— Q'^°°*^7'rootg- 

The constraint C appears in the action with an undetermined multiplier e. The momenta are related to the original 
CT variables by QP^Tphi = Q'°°^^Vroots = (2/3) ^5'^*'^ and T'l^ttric -(l/2)c^^/M'^^ These relations can be used 
to rewrite the Hamiltonian constraint H14|) in terms of the canonical variables. The "smeared" constraint J d^xeC 
generates an infinitesimal conformal transformation Q through the Poisson brackets with ^ = 1 — e/6. 

Let us return to the action as expressed in Eq. 113() . The equations of motion obtained from this action are 
redundant, as shown by Eq. H19|l . because the CT variables are not unique. Said another way, the action (|13(l does 
not depend on gab, Aab, and t separately but only on the combinations 



gab = g ^^^gab = h ^^^hab , 

^ ^ g''^\^h'''\ 

^ab ^ g5/6^a6 ^ _2;il/3 (pah _ p^^ah/g^ ^ 

T = T = P/Vh . 



(21a) 
(21b) 
(21c) 
(21d) 



These quantities are invariant under the conformal transformation |(2J). I will refer to gab, 'f^, A°'^ and f as the 
"invariant CT variables." Note that gab is the unit determinant metric. It is a type (2) tensor density with weight 
—2/3. The variable A°'^ is a traceless type (q) tensor density with weight 5/3. It is proportional to the trace-free 
part of the extrinsic curvature. The variable (p is a, scalar density with weight 1/6, and f is the trace of the extrinsic 
curvature with no density weight. Note that indices on A"'' and other "tilde" quantities are raised and lowered with 
(jab and its inverse 

The invariant CT variables (the variables with tilde's) are combinations of the CT variables (the variables without 
tilde's) that are invariant under the conformal transformation But they can also be viewed as the CT variables 
transformed by Eqs. Q with ^ = g^^/^^. This implies that the action, which is conformally invariant, has the same 
expression in terms of the invariant CT variables as it does in terms of the original CT variables. Thus, we can rewrite 
Eqs. (fU^ . (fTIjl . and (fH^ by placing tilde's on each of the CT variables: 



S^s)[g^^^Aab,^,f,a,p'^] 



dtd-^x 



-]^A''bd^~gab + \fd^g>'' 



(22) 



where 



n 



(9) 



\A^''Aab-l^''r^-l^'R + 4^'D^^, 



Mi^^ = DbA'a - -^'^baf 



(23a) 
(23b) 



Note that the determinant of cjab has been set to one. The full Hamiltonian and momentum constraints arc Ti. = 
H^s'i + and Ma = Mi^^ - ip^ja, respectively. 

In the Hamiltonian density Ti.^^\ the term R is constructed from gab according to the usual formula R = 
2gab ^Qj^f + f Q[bf where t"^^ = g'"^ {dbgdc + dcgbd — ddgbc) /2. Likewise, the Laplacian operator acting on (p 



is defined by D (p — g I dadb^p — r^^9c<^ ) . Because the metric gab itself carries a density weight it transforms under 
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a change of spatial coordinates with an extra factor of J~^/^ as compared with a weight-zero metric. Here, J is 
the Jacobian of the transformation. As a consequence of the density weighting on gab, the terms R and D^ip are 
not separately scalars under spatial coordinate transformations. However, the combination (p~'^R — %(p^^D'^(p, which 
equals the physical scalar curvature 7?., is a scalar. Then together the terms —(p^R/2 + A(p^ D'^ip that appear in the 
Hamiltonian density (|23a|l transform as a scalar density of weight 2. Recall that the lapse anti-density a carries 
weig ht -1, so an^9^ is a weight 1 density as it should be. 

Also observe that the Laplacian acts on (p in Eq. H23a|l as if (p were a scalar rather than a scalar density of 
weight 1/6. This is because the difference between Da acting on a scalar and acting on a scalar density is a term 
proportional to dafj, which vanishes because gab has unit determinant. Likewise, in the momentum density (|23b|) . the 
covariant derivative Db acts on A^a if A\ were a type (]^) tensor with no density weight. 



IV. EQUATIONS OF MOTION FOR THE CONFORMAL-TRACELESS VARIABLES 

Let us derive the dynamical equations of motion by varying the action H22|) with respect to the invariant CT 
variables. Note that the variations and Scjab are traceless. As discussed in the previous section, the functional 
derivatives with respect to A°-^ and gab are trace free: 



5S^a) 1 , s TF 



TF 



Here, iJ^^^ = J dP'xaliS^^ is the gravitational field contribution to the Hamiltonian where HS^^ is given by Eq. (|23a|l . 
The terms d±gab and SH^^^SA"-'' = aAab in Eq. (|24b|l are already traceless, so the equation of motion obtained by 

/ . \TF 

varying the action with respect to is d±gab = —2aAab- The term id^A"^^] can be written as 

TF _ 1 



[d^A-') = d^A-" - -r'-gcdd^A-'' 



O^A'^b _ ^ag-b^cd^^^ ^ (25) 



where the equation d±gab = —'2aAab has been used. The evolution equation for A'^b is found by setting SS^^'^ /Scjab 
equal to —SS'^'"^^ / 5gab = — (a<P^^/2) (s"'') in Eq. H24a|l and using the result from the calculation (|25|l . The remaining 
equations of motion 6S/5(p = and 6S/6f = are straightforward to derive. The complete set is 



d±gab = -2aAab , (26a) 

tTF 

diA"'' =^ 2aA'"=A''-8ip'^ ' 



+ b^^a &'^(p + -(pb°-b^a - -a(pR'''' + -aCp^s"'' 



(26b) 



d^ip = -^aip^T , (26c) 

di_T = ol'^t'^ + a((?R - \\a.(p&(p - ((?&a - X^ipb"- (pb aO. ~ A2ab''ipba<p + ]^a'f{s - 3p) . (26d) 

Note that the matter variables have not been conformally scaled, and the indices on the spatial stress tensor are 
lowered with the physical metric; in particular, s = s°'^hab- 

The equations of motion for the original set of CT variables, gab, Aab, 'P, and t, can be found by cxtremizing the 
action (|13|l . These equations are redundant, as implied by the relation in Eq. (|19() . Alternatively, we can obtain the 
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independent equations of motion by inserting the definitions (I21II into Eqs. (|26|) : 



d±gab 

d±Aab 

d±(p 
d±T 



gffafc d± In, 9 - 2ay^Aab , 
- ^Aab d± In g - 2ay/gAacAl 



aDaDbf + D(a(^ Db)f + ^ipDaDba - ^aipRab + ■^acpsab 



TF 



(27a) 

(27b) 
(27c) 
(27d) 



These equations are essentiahy identical to Eqs. I|26|) . apart from various factors of g and the absence of tilde's. Note 
that the indices are up in Eq. (j26b() . and down in Ec^. l|27b|l . This difference gives rise to the sign difference between 
the first term on the right-hand side of Eq. (|26b|l and the second term on the right-hand side of Eq. I|27b|) . 

It might not be obvious that the terms involving the Ricci tensor Rab and covariant derivative Da on the right-hand 
sides of Eqs. Il26b|) and (|26dp will simplify to the corresponding terms involving Rab and Da on the right-hand sides 
of Eqs. I)27b(l and l|27d|l . However, the following argument shows that this must be the case. The terms involving Rab 
and Da on the right-hand sides of Eqs. (|27|1 can be obtained from the functional derivatives of 



F 



-- I cPxahn 



/ d^xag{(p''R-8ip'D'''tp) 



(28) 



Since F is conformally invariant, we can view it either as a functional of gab and 93 or as a functional of the conformally 
invariant variables g^^^g"'' and g^^^^ip. The functional derivatives with respect to gab and ip are defined by 



SF- 



J 



SF , SF ' 
i — Sgab + -T-Sv 
Sgab Sip 



By splitting SF / Sgab into its trace and trace-free parts, we can rewrite this expression as 

TF 



SF 



-1/3 f 

-g ' gacgbd 



cd 



s{g 



1/12, 



(29) 



(30) 



Then by explicit calculation, we find 



-g ^^^gacgbd 



SF 

Sg, 



cd 



TF 



= [aipRab - ipDaDba - 8D(^aaDb)P - SaDaDbip] 



TF 



5(gi/iV) ^ Sip 



[-a^p'^R + UapD'^ip + A2aD''ipDa'p + lAipD'^ipDaa + <p^D'^a\ 



(31a) 



(31b) 



Since F, g^^^g°'^^ and g^^^^tp are conformally invariant, the expressions on the right-hand sides of Eqs. H31a() and 



(|3rbp must be conformally invariant. Therefore these expressions are unchanged if we conformally transform the 
CT variables to the invariant CT variables. We do this by applying the transformation ^ with ^ = g"^/^^. This 
argument shows that we can place tilde's on the right-hand sides of Eqs. 131|l without changing the values of these 
expressions. Armed with this result, it is straightforward to show that Eqs. H27() follow from Eqs. (|26() . 

The CT equations H27() extremize the action and are equivalent to the ADM equations They are invariant 
under the conformal transformation iPJ. Thus, these equations do not determine the evolution of g, the determinant 
of the conformal metric gab- For practical (numerical) calculations, it can be useful to fix g. As discussed in the 
introduction, a common choice is to specify g = 1. However, we are free to choose g to be any t-dependent spatial 
scalar density of weight 2. There are two natural cases to consider for the evolution of g. The first case is d±g = 0; 
that is, g can be chosen to be constant along the normal to the spacelike hypersurfaces in spacetime. I will refer to 
this as the Eulerian condition, since g is constant for the observers who are at rest in the spacelike hypersurfaces. The 
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second case is dg/dt = 0; that is, g is chosen to be constant along the time flow vector field in spacetime. I will refer 
to this as the Lagrangian condition, since g is constant for the observers who move along the "flow lines" defined by 
the spatial coordinates. Note that the choice g = 1 is a special case of the Lagrangian condition. 
For the Eulerian case the terms d± lug vanish in the equations of motion (|27|l . Thus, we have 

dgab/dt = 2D^,j3h) - 2a^Aat , (32a) 

TF 

(32b) 



aDaDbtp + D(^aCtDh-)ip + ^LpDaDba - ^aipRab + ^afSab 



1 



dgab/dt = 


1 right 


dAab/dt = 


{right 


d(f/dt = 


{right 


dr/dt = 


(right 



d^/dt = /3'=£'c<yJ - -ay/gif'^T , (32c) 

dr/dt ^ P^DcT + aip^r^ + atp^R - Ua(fD'^ip - tp^D'^a - UipD'^ipDaa 

-42aD''ipDaip+^aip^{s-3p) , (32d) 

when the condition d±g — holds. For the Lagrangian case dg/dt = we have 

2 

Qd side of Eq. \ - -gabDcli" , (33a) 

Qd side of Eq. (l32bl) | + \AabDcP'' , (33b) 

Qd side of Eq. (|21cl)| + ^ipD^fi" , (33c) 

Qd side of Eq. (|32d|) | . (33d) 

By contracting with g°-^ we find that Eq. (|32a|l preserves the Eulerian condition. Likewise, Eq. (|33a|) preserves the 
Lagrangian condition. 

The terms 2D(^aPb) in Eqs. H32a() and (|33a|l come from the Lie derivatives of the type (2) tensor gab- The extra term 
— {2/'i)gabDcP'^ that appears in the Lagrangian case H33a|) combines with 2D(^g^/3f,-f to give the same result one would 
obtain by computing the Lie derivative as if gab were a type (2) tensor density with weight —2/3. Similarly, the extra 
terms that appear in the equations of motion for Aab and ip in the Lagrangian case combine with the Lie derivatives 
^/3Aab and Cpip to give the same results one would obtain by computing those Lie derivatives as if Aab were a type 
(2) tensor density of weight 1/3 and ip were a scalar density of weight 1/6. 

As discussed in the introduction, there is no conceptual advantage in leaving the trace of Aab unspecified. If yl = 
were not imposed from the beginning, the independent equations of motion would be identical to Eqs. H27|l with the 

replacements r r + p~^A and Aab (Aab)"^^ ■ The resulting equations would not determine the evolution of r 
and A separately, but only the "trace invariant" combination r + p^^A. We could choose A to be any f-dependent 
spatial scalar. In practice it might be convenient to determine A by removing the trace-free (TF) symbol from the 
terms in square brackets on the right-hand side of Eq. (|27bp . 



V. BSSN EQUATIONS 

The conformally invariant BSSN system of evolution equations 0,0 is based on the CT equations (|27|l . It uses 
new fields, the "conformal connection functions" , defined by 

T^^g'^'^ri, = ~l^da{V99n ■ (34) 

Here, FJ^j, are the Christoffel symbols built from the conformal metric gab and da stands for the partial derivative with 
respect to the spatial coordinate In the original treatment of the BSSN system, the condition g = 1 was imposed 
from the outset. As in the previous sections, I do not impose g ~ I. 

The key idea of the BSSN system is to replace certain derivatives of the metric that appear in the Ricci tensor on 
the right-hand side of Eq. Ij27b|l by the conformal connection functions F". By explicit calculation, we have 

Rab = --^g'^'^dcddgab + gciadb)^" + r''F(ab)c + g""^ (^^^l{a^b)ed + ^IJ" ebdj , (35) 
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where Tabc = gadXfc- The CT system (|27|) must be extended to include an evolution equation for F". This is found 
by computing d^T"" from its definition H34|l . For this calculation, it is helpful to introduce the invariant conformal 
connection functions, F"^ = .9"''F^^ = —dafj'^'^- These variables are invariant under the conformal transformation (j^J 
since they are built from gab = g^^^'^gab- They are related to F" by 

pa ^ ^l/Spa ^ ^l/Q)g-V3g-b9^g . (36) 

The evolution equation for F° is straightforward to compute once we recognize that d± and da commute.^ The result 
can be expressed as 

d± [g'/'r-) = (gi/^^^^^lns) - 24 (^ag'^'^A-') , (37) 

where dblng — g^^dbg. 

The equations H27() . H37() appear different from the familiar expression of the BSSN equations largely due to differ- 
ences in notation. Let me define a set of "BSSN variables" , denoted by carets: 

ip = e"^ , (38a) 

Aab = e^^Aab , (38b) 
a EE e-^^a/^ . (38c) 

The definitions for the metric gab, trace of extrinsic curvature r. and shift vector /3° are unchanged. Under a conformal 
transformation, the BSSN variables change by 

(39a) 
(39b) 
(39c) 
(39d) 
(39e) 

The scalar lapse function a and the shift vector /3° are conformally invariant. 
With the change of notation the BSSN system l(?7|l . (p!7jl becomes 

d^gab = ^gabd^lng -2aAab , (40a) 
d±Aab = -^Aabd^\ng -2aAacAl + ^AabT 

+e-^^ [-2aDaDb(p + AaDa^Dbip + ADf^aaDb) (p - DaDbOt + aRab ~ asab] '^^ , (40b) 

dy_(p = --^ 9_L In 5 - ^ar , (40c) 

12 D 

d^T = dr^ + e~^^ [aR - SaD^p - SaD^pDaP - D'^a - 2D°'aDaP>) + ^a(s - 3p) , (40d) 
diT'^ = -^T'' d^lng - ^g'^^'dbd^lY^g ^ ^db(a^A'^'') . (40e) 

We can further modify these equations by making use of the Hamiltonian and momentum constraints HS^^ = —gif^^p 
and A^i^"* = \/9^^ja^ where Ti.^^^ and A^i^' are given by Eqs. p4|) and p8|) . respectively. First, we rewrite the 
momentum constraint as 

db [Vae'^A^') = ^e'^ (^-nJ"' + ^g'^'dbT + g^'jb^ , (41) 



9ab - 


K gab , 


(p ~ 


-> ^-In^ , 


Aab - 


- £,^Aab , 


T — 


T , 


pa 





That is, £^ and 9a commute. To show this, recall that the components of, say, a contravariant vector in a "barred" coordinate system 
a;" are related to the components of that vector in an unbarred coordinate system by v'^{x) = {x)di,x°' . The coordinate derivatives 
of the vector components are related by daV^{x) = {dax'^}dclv'^ {x)ddS^]- Similar expressions hold for the components, and derivatives 
of components, of other types of tensors. Define the Lie derivative by CjsT = T{x) — T{x) where = — /3° and f) is infinitesimal. 
(Indices have been suppressed on T.) One then finds that Cf3{daT) = da{CgT). 
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and use this result to replace the spatial derivatives of A"'^ on the right-hand side of Eq. (|40e(l . Next, we rewrite the 
Hamiltonian constraint as 



R 



r 



SiD^^ + D'^ifDaip) 



(42) 



and use this result to replace the Ricci scalar R on the right-hand side of Eq. Ij40d|l . The end result of these changes 
is 



d±gab 

d±Aab 

d±ip 
d±T 



1 



gab d± In g - 2a Aab , 
-^iat d± Ing - 2aAacAl + aAabT 

-e'^v [^2aDaDbip + AaDacpDbip + ADi^aaDt^ip - DaDbCt + aRab 



asabl 



iTF 



12 



d± In g 



-ar 



aAabA-' 



e-4v (•£)2^ ^ 2D''aDa0) + -a{s + p) 



-T" di_ In.g - -g''''dbd^ Ing- 2A''^dba + 2a 
3 ■ 6' 



eA-'dbip + TlA'"^ - -g-'dbT - g-'jt 



Equations H43|l with the Ricci tensor Rab given by Eq. (|35|) define the conformally invariant BSSN system. 
We can rewrite the BSSN equations in a more compact form by using the identities 

V^d = e"'''^ {D^a + 2D''aDaip) 



in Eq. (j43d|l and 



aTZab - VaVbd] = [-2aDaDbip + AaDa(pDb(p + AD^aC^Di^-^Lp - DaDba + aRab] 



TF 



(43a) 

(43b) 
(43c) 

(43d) 

(43e) 

(44) 
(45) 



in Eq. (|43b|l . Recall that TZab and Va arc the Ricci tensor and covariant derivative constructed from the physical 
metric hab = e'^'f'gab- 

The BSSN equations (|43|l are invariant under the conformal transformation H39|l . As with the CT equations H27|l . 
these equations do not determine the evolution of g. We must specify d±g, which appears in several places on the 
right-hand sides of Eqs. H43|l . as a separate condition. Alternatively, we could write the BSSN system in terms of 
conformally invariant variables, including the invariant conformal connection functions F". This would remove the 
terms proportional to d± In g in the BSSN equations but add certain density weights to each of the variables (except 
r). As defined here, gab and Aab are type (2) tensors with no density weight, and (f and r are scalars with no density 
weight. The coordinate transformation rule for the conformal connection functions follow from their definition H34() 
and the familiar inhomogeneous transformation rule for the Christoffel symbols. In particular, F" carries no density 
weight. 

Two natural choices for the evolution of g are the Eulerian condition d±g = and the Lagrangian condition 
dg/dt = 0. With the Eulerian condition, we have 



dgab/dt 


= '^D(aPb) - 


- 2aAab , 






(46a) 


dAab/dt 


= P^D.Aab 


+ 2i,(,Z?b)/3^ - 2aAacAl 


+ aAabT + 


s ^"^ [aTZab - '^a'^ba - aSab]^^ , 


(46b) 


d'f/dt 




1 






(46c) 


dr/dt 


= I3'D,T + 


idr^ + diabi"' - V^d 4 


1 

-a{s + p) 




(46d) 


dV'/dt 


= P'^dcT" - 


r'^dcP'' + g'^'dbdcP" - 2A 


^^dba + 2d 


QA-^db^ + Vt^A"- - |g'"'a,T - g'^^'jb 


. (46e) 



The identities and have been used to express these equations in compact form. For the Lagrangian condition. 
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we have 



dgab/dt = 


I right-hand side of Eq. H46a|) > — 


dAab/dt = 


< right-hand side of Eq. I)46b|) > — 


dip/dt = 


< right-hand side of Eq. H46c|l > + 


dr/dt = 


< right-hand side of Eq. Ij46d() > , 


dV'/dt = 


< right-hand side of Eq. H4t)e|l > + 



(47a) 
(47b) 
(47c) 
(47d) 
(47e) 



The Eulerian and Lagrangian conditions are preserved by equations 146a|) and (|47a|l . respectively. Note that the 
restriction g = 1 is a special case of the Lagrangian condition. Also note that with g = 1, the extra terms in Eqs. H47(l 
simplify since in that case Dc(3'' = dcP'^- 

It has been suggested that the "Gamma freezing" condition, dT°- /dt = might be useful as a means of speci fying 
the shift vector for numerical evolutions based on the BSSN system. The related "Gamma driver" conditions |17l |. 
which are less time-consuming to solve numerically than the Gamma freezing condition, have been used with some 
success. Here we note that these conditions depend specifically on the way that the shift terms enter the evolution 
equation for F". If we choose the Eulerian condition d±g = to break the conformal invariance, then the Gamma 
freezing shift equation is 



g'"^dbdc(i°' = terms containing at most first derivatives of . 



(48) 



If we choose the Lagrangian condition dg/dt = to break conformal invariance, the Gamma freezing condition 
becomes 



1 



g dbdcP"' + -g dbdcP'^ = terms containing at most first derivatives of /3° . 



(49) 



It might be interesting to explore the difference between these two shift conditions. 



VI. CONFORMAL THIN SANDWICH EQUATIONS 

Initial data f or g eneral relativity must satisfy the constraint equations. In the original York-Lichnerowicz conformal 
decomposition [l^ , the gravitational field parts of the constraints are written as in Eqs. (|14|l and (|18() . The Hamiltonian 
constraint -ip^'^gp can be solved for the Laplacian of the conformal factor (p, 

D'v = -\ip~'A^'Aa, + + l^i? - . (50) 

The trace free extrinsic curvature Aab is split into a transverse part and a longitudinal part. With the so-called 
"conformal transverse-traceless decomposition" , the longitudinal part of Aab is expressed in terms of derivatives of a 
vector X". Then the principal part of the momentum constraint Ma = '-P^ \fgja is proportional to 

AtXa = Db(lLX)\ , (51) 

where 

{-!LX)ab = 2i?(,Xb) - (2/3).gahi5,X= . (52) 

Ai is an elliptic operator. In this way, the Hamiltonian and momentum constraints are expressed as a system of 
elliptic equations for and The freely specifiable parts of the gravitational field are the conformal metric gab^ 
the trace of the extrinsic curvature r, and the transverse part of Aab- 

More recently, York |^ has recognized that the momentum constraint can be expressed in terms of the ellip- 
tic operator acting on the shift vector. Observe that the trace-free part of the conformal metric velocity is 
{dgab/dt)^^ = dgab/dt — {gab/'i) d{lng)/dt. Using the "dot" notation for time derivatives, we find that Eq. H27a|) 
becomes 

Aab = -^{gl^~imab) (53) 
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where {Ej(3)ab is defined in Eq. H52() . Inserting this result into the momentum constraint, we find 

Ai/3a = D'-gJ^ + 2^AabD''a + ^a^^^D^r + 2a^^'^ja , (54) 

where is defined by Eq. H51|l . As described in Ref. 0, one can specify freely the gravitational quantities gab, iiab , 
T, and a then solve Eqs. H5U|) and (|54|) for ip and Wherever A"''' appears, it is written in terms of the conformal 
metric and its derivatives via Eq. (jSHfl . 

The conformal thin sandwich construction Q is an extension of this analysis to include the lapse anti-density as one 
of the unknowns. By using the Hamiltonian constraint (|50|1 to eliminate the Laplacian of (p, we find that Eq. (j27d|l 
can be written as 

D^a = -—^-^d±T - -a^^T^ - -aR + -a^-"" A''''Aab 
Va 6 4 4 

-Uip-^D^'ipDaa - A2aip-^D''ipDa(p + ^aip'^is + 4p) . (55) 

The principal part of this equation is the conformal Laplacian operator acting on a. For the conformal thin sandwich 
initial data construction the freely specified quantities are the conformal metric gab, the trace-free part of the metric 
velocity c/ab, the trace of the extrinsic curvature r, and its time derivative f = dr/dt. Equations (|50|1 . (|54|l . and 
(|55|l . along with the expression for A"''^ given in Eq. H53|l . constitute an elliptic system of equations to be solved for 
the conformal factor the shift vector and the lapse anti-density a. Once these equations are solved, the initial 
data in terms of the physical metric hab and extrinsic curvature Kab can be obtained from Eqs. 

The conformal thin sandwich equations H5U|) . H53|) . (|54|) . and (|55|) are derived from the CT equations and the 
Hamiltonian and momentum constraints, all of which are conformally invariant. It follows that the conformal thin 
sandwich equations are invariant under the conformal transformation |(5J. To be precise, let and a denote the 

solution of the conformal thin sandwich equations for a given set of input data gab, sJfT, and f. Then the solution 
for input data ^'^gab, ^^ffJh^: t, and f will be (fi/^, /3", and a. The physical metric hab and extrinsic curvature Kab are 
the same in either case. 

In evolving the conformal thin sandwich data one can choose the lapse anti-density and shift vector a and freely, 
without regard to the values obtained from the initial data construction. However, if the initial values of the lapse 
and shift are chosen to coincide with the values of a and /3° that were computed from the thin sandwich equations, 
then initially the trace-free part of the conformal metric velocity will coincide with the chosen value of gj^ and the 
initial time derivative of the extrinsic curvature's trace will be given by the chosen value of r. This will be the case 
whether the data is evolved via the ADM equations ||SJ), the CT equations (|27|) . or the BSSN equations H43|l . In the 
later two cases, we are free to choose the Eulerian or Lagrangian condition, or any other condition for the evolution 
of g. 



VII. CFC EQUATIONS 

The approximate evolution equations obtained from the CFC have been used in numerical studies of binary neutron 
star coalescence @, and supernovae 0, |23| . This approximation to general relativity appears to be quite good , 
at least for systems that are not too far from spherical symmetry. 

It turns out that the CFC equations are precisely the conformal thin sandwich equations, lf3n|) . and 
The CFC approximation can be described in a generalized way as follows. Let the conformal metric gab and the trace 
of the extrinsic curvature r be specified freely for all times. Then gab, gJb ■• '''i ^^'^ known for all times and 

the thin sandwich equations can be solved for the gravitational data a, /3°, if, and Aab as soon as the source data p, 
ja, and s"'' is available. The idea, then, is to solve the conformal thin sandwich equations at the initial time using 
the initial values for the sources. We then evolve the sources forward in time to the first timestep using the matter 
equations of motion. These equations depend on the physical metric hab ~ '~P^9ab, lapse anti-density, and shift. From 
the source values at the first timestep, the conformal thin sandwich equations are solved to complete the gravitational 
data at the first timestep. The process repeats. 

For the CFC approximation, one sets the conformal metric equal to a flat metric for all times. In particular g^ ~ 
in Eq. H54|) . Also, maximal slicing is assumed so that r = for all times. Recall that the conformal thin sandwich 
equations are equivalent to the Hamiltonian constraint, the momentum constraint, and two of the four CT evolution 
equations, namely, Eqs. (|27a|l and Ij27d|l . Thus, the approximation in the CFC formalism consists in ignoring the 
evolution equations for Aab and (p, namely, Eqs. Ij27b|) and H27c|) . 
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VIII. SUMMARY 



The conformal-traceless decomposition of the gravitational field appears in a number of contexts in general relativity, 
most notably in the analysis of the initial value problem and in the construction of the BSSN system of evolution 
equations. In this paper I have written the action functional and equations of motion in terms of conformal-traceless 
variables. I do not invoke the frequently imposed condition g = 1 on the determinant of the conformal metric. As 
a consequence, the action and equations of motion are conformally invariant. I have presented two possibilities for 
breaking conformal invariance, namely, the Eulerian condition d^g = and Lagrangian condition dg/dt = on 
the evolution of g. I also extended the equations of motion to obtain a conformally invariant version of the BSSN 
system. For this system as well, the invariance can be broken by specifying the Eulerian or Lagrangian condition. I 
showed that the conformal thin sandwich equations for gravitational initial data are obtained from the Hamiltonian 
and momentum constraints along with a subset of the evolution equations written in conformal-traceless variables. 
Finally, I have pointed out that the CFC approximation to general relativity consists in solving the conformal thin 
sandwich equations at each time step, assuming the conformal metric is flat and the time slicing is maximal. 
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